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Yale University, New Haven, Connecticut 06511, USA 

We present measurements of the time-dependent CP-violation parameters S and C in rj'K'^ 
decays. The data sample corresponds to 384 million BB pairs produced by e+e~ annihilation at 
the T(45). The results are S = 0.58 ± 0.10 ± 0.03, and C = -0.16 ± 0.07 ± 0.03. We observe 
mixing-induced CP violation with a significance of 5.5 standard deviations in this b s penguin 
dominated mode. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Measurements of time-dependent CP asymmetries in 
B° meson decays through Cabibbo-Kobayashi-Maskawa 
(CKM) favored b ccs amplitudes [H have provided cru- 
cial tests of the mechanism of CP violation in the Stan- 
dard Model (SM) 0. Decays of B° mesons to charmless 
hadronic final states such as rj'K'^ proceed mostly via a 
single loop (penguin) amplitude. In the SM the penguin 
amplitude has approximately the same weak phase as 
the b CCS transition, but it is sensitive to the possible 
presence of new heavy particles in the loop Q . The mea- 
surement of CP asymmetries in —> j]' thus provides 
an important test for such effects. 

Within the SM, CKM-suppressed amplitudes and mul- 
tiple particles in the loop introduce additional weak 
phases whose contribution may not be negligible [1, @, 0, 
|7|. The time-dependent CP- violation parameter S (de- 
fined in Eq.[T] below) measured in the decay -q' 
is compared with the value of sin2/3 from measurements 
of time-dependent CP violation in B decays to states 
containing charmonium and a neutral kaon. The devia- 
tion AS = S — sin2/3 has been estimated in several the- 
oretical approaches: QCD factorization (QCDF) Q, 
QCDF with modeled rescattering Soft Colhnear Ef- 
fective Theory 0, and SU(3) symmetry These 
models estimate |AS'| to be of the order 0.01, and with 
uncertainties give bounds \A.S\ ^ 0.05. 

The time-dependent CP asymmetry in the decay 
B^ v 'K'^ has been measured previously by the 
BaBAR and Belle [3 Collaborations. In this Let- 
ter we update our previous measurements using an inte- 
grated luminosity of 349 fb~^, corresponding to 384 ± 4 
million BB pairs, recorded at the T{AS) resonance 
(center-of-mass energy_ys = 10.58 GeV). Belle has since 
updated their resul ts [ l4l| . Our data were collected with 
the BaBAR detector [15| at the PEP-II asymmetric-energy 
e+e^ colhder. In addition to the B° — » r]'K'^ decays used 
previously, we now also include the decay i?" v'^^l- 

Charged particles from e'^e~ interactions are detected, 
and their momenta measured, by a combination of five 
layers of double-sided silicon microstrip detectors and 
a 40-layer drift chamber, both operating in the 1.5 T 
magnetic field of a superconducting solenoid. Photons 
and electrons are identified with a CsI(Tl) electromag- 
netic calorimeter (EMC). Charged particle identification 
is provided by the average energy loss in the track- 
ing devices and by an internally reflecting ring imaging 
Cherenkov detector covering the central region. The in- 



strumented flux return (IFR) of the magnet allows the 
identification of muons and mesons. 

We reconstruct a decaying into the CP eigenstate 
Ty'ifg or rj'K^ (Bcp)- From the remaining particles in the 
event we also reconstruct the decay vertex of the other 
B meson (Stag) and identify its flavor. The difference 
At = tcp — ttag of the proper decay times tcp and itag 
of the CP and tag B mesons, respectively, is obtained 
from the measured distance between the Bcp and Stag 
decay vertices and from the boost (/?7 = 0.56) of the 
e'^e~ system. The At distribution is given by: 



F(At) 



HAt|/r 



4t 



[1 T Aw ± 



(1) 



(1 - 2w) {-j]Ssm{AmdAt) - C cos(AmrfAi))] 



where rj is the CP eigenvalue of the final state (—1 for 
rj'K]^, +1 for rj'K^). The upper (lower) sign denotes a 
decay accompanied by a 5° (S°) tag, r is the mean B° 
lifetime, Ato^ is the mixing frequency, and the mistag 
parameters w and Aw are the average and difference, 
respectively, of the probabilities that a true is incor- 
rectly tagged as a i?'^ or vice versa. The tagging algo- 
rithm has six mutually exclusive tagging categories and 
a measured analyzing power of (30.4±0.3)% [16|. A non- 
zero value of the parameter C would indicate direct CP 
violation. 

We establish the event selection criteria with the aid of 
a detailed Monte Carlo (MC) simulation of the B produc- 
tion and decay sequences, and of the detector response 



17| . These criteria are designed to retain signal events 



with high efficiency while removing most of the back- 
ground. 

The _B-daughter candidates are reconstructed through 
their decays tt" -> 77, ?? ^ 77 iV'yf)^ V ^ 7r+7r~7r° {r]^^), 
r,' Ty^-yTT+TT- {ii'^f^^^^^J, 77' 

/ n / / \ 

n 



/0°7 iVp-y), where p° 



+ 7r-, K° TT+TT- 

(if°+,-) or ^-0^0 (i^Oo^o). Only the r^'^^^^^^^ mode is 
used for the rj'Kl sample. The requirements on the 
invariant masses of these particle combinations are the 
same as in our previous analysis [HI . The list of all decay 
modes used in the current analysis can be seen in Table 
m Signal candidates are reconstructed from clusters 
of energy deposited in the EMC or from hits in the IFR 
not associated with any charged track in the event [l^]- 
From the cluster centroid and the i?° decay vertex we 
determine the direction (but not the magnitude) of the 
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momentum p^o . 

For ry'ifg decays we reconstruct the B-meson candidate 
by combining the four-momenta of the ifg and 77' with 
a vertex constraint. We also constrain the rj, rj' , and tt*^ 
masses to world-average values jl9l] . From the kinemat- 
ics of T{4S) decays we determine the energy-substituted 

mass Toes = \J (^s -I- po • Pb)'^/Eq — and the energy 

difference AE = Eg — ^^/s, where (i?o,Po) and {Eb,Pb) 
are the laboratory four- momenta of the T{AS) and the 
B candidate, respectively, and the asterisk denotes the 
T(4S') rest frame. The resolution is 3 MeV in toes and 
20 — 50 MeV in AE, depending on the decay mode. 

For rj'Kl candidates we obtain AE and p^o from a 
fit with the and masses constrained to world- 
average values |19| . To make a match with the measured 
K'l direction we construct the missing momentum Pmiss 
from Po and all charged tracks and neutral clusters other 
than the candidate. We then project Pmiss onto p^o , 
and require the component perpendicular to the beam 
line, Crii^ to satisfy pfJi^ - p^o^ > -0.5 GeV.^ This 
value was chosen to minimize the yield uncertainty in the 
presence of background. 

For r/K^ we require 5.25 < toes < 5.29 GeV and 
\AE\ < 0.2 GeV, for r]'K° we require -0.01 < AE < 
0.04 GeV, and for aU decays \At\ < 20 ps, and, for the 
error on At, a^t < 2.5 ps. 

Background events arise primarily from random com- 
binations of particles in continuum e^e^ —^ qq events 
{q — u,d,s,c). We reduce these with requirements on 
the angle Ot between the thrust axis of the B candidate 
in the T(4S') frame and that of the rest of the charged 
tracks and neutral calorimeter clusters in the event. In 
the fit we discriminate further against qq background 
with a Fisher discriminant J-^ that combines several vari- 
ables that characterize the production dynamics and en- 
ergy flow in the event (20l | . For the 77^^ decays we require 
I '^os^dccl ^ t° reduce the combinatorial background. 
Here 6^^^ is the angle between the momenta of the 
daughter tt" and of the 77', measured in the rest frame. 

For B'-' r}' candidates we require that the cosine 
of the polar angle of the total missing momentum in 
the laboratory system be less than 0.95, to reject very 
forward qq jets. The purity of the K^^ candidates re- 
constructed in the EMC is further improved by a re- 
quirement on the output of a neural network (NN) that 
takes cluster-shape variables as inputs. The NN was 
trained on MC signal events and data events in the re- 
gion 0.02 < AE < 0.04 GeV. We check the performance 
of the NN on data with candidates in the larger 
B° J/ipKl data sample. 

The average number of candidates found per selected 
event is between 1.08 and 1.32, depending on the final 
state. In the case of events with multiple candidates we 
choose the candidate with the smallest value of a con- 
structed from the deviations from expected values of one 



or more of the daughter resonance masses, or with the 
best decay vertex probability for the B, depending on 
the decay channel. Furthermore, in the rj' sample, if 
several B candidates have the same vertex probability, 
we choose the candidate with the information taken 
from, in order, EMC and IFR, EMC only, or IFR only. 
From the simulation we find that this algorithm selects 
the correct-combination candidate in about two thirds of 
the events containing multiple candidates. 

We obtain the common CP-violation parameters and 
signal yields for each channel from a maximum likelihood 
fit with the input observables AE', toes, and At. The 
selected sample sizes are given in the first column of Ta- 
ble m We estimate from the simulation a contribution to 
the input sample of less than 1.1 % of background from 
other charmless B decay modes. These events have final 
states different from the signal, but similar kinematics, 
and exhibit broad peaks in the signal regions of some ob- 
servables. We find that the BB background component 
is needed only for the channels with -q'p^. We account 
for these with a separate component in the probability 
density function (PDF). For each component j (signal, 
qq combinatorial background, or BB background) and 
tagging category c, we define a total probability density 
function for event i as: 

Kc ^ V,{m^sl-V,{AE^)-T,{n-V,{At\ a\,;c) , (2) 

except for rj' K^^ for which ^^(toes*) is omitted. The fac- 
tored form of the PDF is a good approximation since 
linear correlations are small. 

We write the extended likelihood function for all events 
of the decay mode d as 



(3) 



where Uj is the yield of events of component j, fj^c is 
the fraction of events of component j for each category 
c, ric = nsigfsig,c + riqqfqq.c + nBBfBB,c is the number 
of events found by the fitter for category c, and Nc is 
the number of events of category c in the sample. When 
combining decay modes we form the grand likelihood C = 
Y\Cd- We fix both fsig^c and fBB,c to /bh.^.c, the values 
measured with the large sample of fully reconstructed 5° 
decays into flavor eigenstates (Bflav sample) [isj . 

The PDF 7'sig(Ai, aAt',c), for each category c, is the 
convolution of F{At; c) (Eq.[T]) with the signal resolution 
function (sum of three Gaussians) determined from the 
-Bflav sample. The other PDF forms are: the sum of two 
Gaussians for Psig('nEs) and Vsig{AE); the sum of three 
Gaussians for Vqq (At; c) and Vgg{At; c); an asym- 
metric Gaussian with different widths below and above 
the peak for Vj{J-) (a small "tail" Gaussian is added 
for Vqq {J-)); a finear dependence for Vqq {AE) and a 
fourth-order polynomial for Vgg{AE); for Vqq (toes) 
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exp [—£(1 — x"^) 



and V gg{m-E,s) the function xy/Y 

with X = 2mEs/\/s and £ a free parameter 2lj and the 
same function plus a Gaussian, respectively. 

For the signal and BB background components we de- 
termine the PDF parameters from simulation. We study 
large control samples of B decays to charm final states 
of similar topology to verify the simulated resolutions in 
AE and mss, adjusting the PDFs to account for any dif- 
ferences found. The qq background parameters are free 
to vary in the final fit. Thus, for the six channels listed in 
Table m we perform a single fit with 93 free parameters: 
S, C, signal yields (6), rj'p^K^ BB background yields 
(2), continuum background yields (6) and fractions (30), 
background At, toes, AE, T PDF parameters (47). The 
par ameters t and Am^ are fixed to world-average values 

We test and calibrate the fitting procedure by apply- 
ing it to ensembles of simulated experiments with qq 
events drawn from the PDF into which we have embed- 
ded the expected number of signal and BB background 
events randomly extracted from the fully simulated MC 
samples. We find negligible bias for C. For S we find 
and apply multiplicative correction factors for bias from 
dilution due to cross-feed from BB background to sig- 
nal events equal to f.03 in the final states ri'K^. 



r;(77)7i 



^Kl and Vp-,Ko.^ 



TABLE I; Results of the fits. Subscripts for r;' decay modes 
denote Jy^f^.^)^^ (1), v'p-, (2), and v'„(an)7V7v (3)- 



Mode 


# events Signal yield 


S 


C 




664 


224 ± 16 


0.61 ±0.23 


-0.26 ±0.14 




11943 


566 ± 30 


0.56 ±0.14 


-0.24 ±0.10 




177 


73 ±9 


0.89 ± 0.35 


0.14 ±0.25 




490 


52 ±9 


0.84 ±0.42 


-0.26 ±0.36 




13915 


133 ± 24 


0.56 ± 0.41 


0.15 ±0.27 








0.62 ± 0.11 


-0.18 ±0.07 




4199 


204 ± 24 


0.32 ± 0.28 


0.08 ±0.23 


r,'K' 






0.58 ±0.10 


-0.16 ±0.07 
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FIG. 1: Distributions projected (see text) onto (a) rriES and 
(b) AE for ri' Ks candidates, and (c) AE for tj'KI candidates. 
The solid lines shows the full fit result and the dashed lines 
show the background contributions. 

Results from the fit for the signal yields and the CP 
parameters S and C are presented in Table [H In Fig. [1] 



we show the projections onto TOes and AE for a subset 
of the data for which the ratio between the likelihood 
of signal events and the sum of likelihoods of signal and 
background events (computed without the variable plot- 
ted) exceeds a mode-dependent threshold that optimizes 
the sensitivity. In Fig. [5] we give the At and asymme- 
try projections of the events selected as for Fig. [1] We 
measure a correlation of 3.2% between S and C in the 
fit. 

We perform several crosschecks of our analysis tech- 
nique including time-dependent fits for B"*" decays 
to the charged final states r]l^^^^^^^K+ , Vp-yK"^ , and 
^r;(37r)ir;r^^' removing One fit variable at a time; fits 
without BB PDFs; fits with multiple BB components; 
fits allowing for non-zero CP information in BB events; 
fits with C ~ 0. In all cases, we find results consistent 
with expectation. The value S = 0.62±0.11 for -q'K'^ dif- 



fers from our previous measurement S — 0.30 ±0.14 [12| 
due to the improved event reconstruction (with a contri- 
bution of +0.08) and selection (±0.12) and to the addi- 
tional data collected (±0.12). With a model of the data 
sample changes introduced by our revised event recon- 
struction and new data, we find that our current result 
has a statistical probability of 35% (50%) for an assumed 
true value of S of 0.61 (0.70). 
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At (ps) 




FIG. 2: Projections (see text) onto At for (a) r]' and (c) 
rj' K1 of the data (points with error bars for i?" tags (Ngo) 
in red empty rectangles and B" tags {Ngo) in in blue solid 
circles), fit function (red dashed and blue solid lines for 
and tj" tagged events, respectively), and background function 
(black shaded regions) . We show the raw asymmetry, (A^^o — 
Ngo)/{NBO + Ngo), for (b) r/'Jfg and (d) r;'A"2; the lines 
represent the fit functions. 

We have studied the systematic uncertainties arising 
from several sources (in decreasing order of magnitude): 
variation of the signal PDF shape parameters within their 
errors, modeling of the signal At distribution, use of At 
signal parameters from the Sflav sample, interference be- 
tween the CKM-suppressed b — > ucd amplitude and the 
favored b cud amplitude for some tag-side B decays 
[i^ . BB background, SVT alignment, and position and 
size of the beam spot. The i?flav sample is used to deter- 
mine the errors associated with the signal At resolutions, 
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tagging efRciencies, and mistag rates. We take the un- 
certainties in tb and Am^ from published measurements 
[lit . Summing all systematic errors in quadrature, we 
obtain 0.03 for S and 0.03 for C. 

In conclusion, we have used a sample containing 
1252 ±50 flavor-tagged rj' events to measure the time- 
dependent CP violation parameters, S — 0.58 ± 0.10 ± 
0.03 and C = -0.16 ± 0.07 ± 0.03. This sample is 2.1 
times as large as that of our previous measurement [l^ . 
Our result for S is consistent with the world average 
of sin2/3 measured in — > J/tpK^. We observe 

mixing- induced CP violation in decays to -q' with 
a significance (systematic uncertainties included) of 5.5 
standard deviations. Our result for direct-CP violation 
is 2.1 standard deviations from zero. 
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